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Abstract. We present the first results showing that the Frechet distance 
between non-flat surfaces can be approximated within a constant factor 
in polynomial time. Computing the Frechet distance for surfaces is a sur- 
prisingly hard problem. It is not known whether it is computable, it has 
been shown to be NP-hard, and the only known algorithm computes the 
Frechet distance for flat surfaces (Buchin et al.). We adapt this algorithm 
to create one for computing the Frechet distance for a class of surfaces 
which we call folded polygons. Unfortunately, if extended directly the 
original algorithm no longer guarantees that a homeomorphism exists 
between the surfaces. We present three different methods to address this 
problem. The first of which is a fixed-parameter tractable algorithm. 
The second is a polynomial-time approximation algorithm which ap- 
proximates the optimum mapping. Finally, we present a restricted class 
of folded polygons for which we can compute the Frechet distance in 
polynomial time. 

Keywords: Computational Geometry, Shape Matching, Frechet Distance 
1 Introduction 

The Frechet distance is a similarity metric for continuous shapes such as curves 
and surfaces. In the case of computing it between two (directed open) curves 
there is an intuitive explanation of the Frechet distance. Suppose a man walks 
along one curve, a dog walks along the other, and they are connected by a leash. 
They can vary their relative speeds but cannot move backwards. Such a walk 
pairs every point on one curve to one and only one point on the other curve (i.e., 
creates homeomorphism between surfaces) in a continuous way. The Frechet 
distance of the curves is the minimum leash length required for the man and 
dog to walk along these curves. Although less intuitive, the idea is similar for 
surfaces. 

While the Frechet distance between polygonal curves can be computed in 
polynomial time |Alt95] . computing it between surfaces is much harder. In 
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[God98 it was shown that even computing the Frechet distance between a trian- 
gle and a self-intersecting surface is NP-hard. This result was extended in [BuclOj 
to show that computing the Frechet distance between 2d terrains as well as be- 
tween polygons with holes is also NP-hard. Furthermore, while in All l()j it was 
shown to be upper semi-computable, it remains an open question whether the 
Frechet distance between general surfaces is even computable. 

On the other hand, in }Buc06| a polynomial time algorithm is given for com- 
puting the Frechet distance between two (flat) simple polygons. This was the 
first paper to give any algorithm for computing the Frechet distance for a non- 
trivial class of surfaces and remains the only known approach for computing it. 
The main idea of their algorithm is to restrict the kinds of different mappings 
that need to be considered. Our contribution is to generalize their algorithm to 
a class of non-flat surfaces we refer to as folded polygons. Given that theirs is 
the only known approach it is of particular importance to explore extending it 
to new classes of surfaces. The major problem we encountered in generalizing 
the work of [Buc06 was that the kinds of different of mappings which need to 
be considered is less restricted. We address three different methods to resolve 
this problem. In Section [31 we outline a fixed-parameter tractable algorithm. In 
Section 21 we describe a polynomial-time approximation algorithm to compute 
the Frechet distance between folded polygons within a constant factor. In Sec- 
tion [5l we describe a nontrivial class of folded polygons for which the original 
algorithm presented in IBuc06l will compute an exact result. 

2 Preliminaries 

The Frechet distance is defined for two k-dimensional hypersurfaces P,Q : [0,l] x — > 
R d , where x < d, as 

6 F (P, Q) = inf sup \\P(p) - Q(a(p))\\ 

<t:A^B peA 

where a ranges over orientation-preserving homeomorphisms that map each 
point p G P to an image point q = <r(p) G Q. Lastly, || ■ || is the Euclidean 
norm but other metrics could be used instead. 

Let P, Q : [0, l] 2 — > M. d be connected polyhedral surfaces for each of which 
we have a convex subdivision. We assume that the dual graphs of the convex 
subdivisions are acyclic, which means that the subdivisions do not have any 
interior vertices. We will refer to surfaces of this type as folded polygons. We 
refer to the interior convex subdivision edges of P and Q as diagonals and edges 
respectively. Let m and n be the complexities of P and Q respectively. Let k 
and I be the number of diagonals and edges respectively. Assume without loss 
of generality the number of diagonals is smaller than the number of edges. Let 
T m atrixmuit(N) denote the time to multiply two N x N matrices. 

2.1 Simple Polygons Algorithm Summary 

In previous work Buchin et al. |Buc06 compute the Frechet distance between 
simple polygons P, Q. The authors show that, while the Frechet distance be- 
tween a convex polygon and a simple polygon is the Frechet distance of their 



boundaries, this is not the case for two simple polygons. The idea of their algo- 
rithm is to find a convex subdivision of P and map each of the convex regions 
of it continuously to distinct parts of Q such that taken together the images 
account for all of Q. 

First, they show that the decision problem Sp(P,Q) < e can be solved by 
(1) mapping the boundary of P, which we denote by dP, onto the boundary of 
Q, which we denote by dQ, such that Sp(dP,dQ) < e and (2) mapping each 
diagonal d in the convex subdivision of P to a shortest path / C Q such that 
both endpoints of / lie on dQ and such that Sp(d, f) < e. 

In order to solve subproblem (1) they use the notion of a free space diagram. 
For open curves /, g : [0, 1] — > R d it is defined as FS e (f, g) = {(x, y) \ x £ f, y G 
9, \\x — y\\ < e} where e > 0. A monotone path starting at the bottom left 
corner of the free space diagram going to the top right exists if and only if the 
curves are within Frechet distance e. As shown in |Alt95j . this can be extended 
to closed curves by concatenating two copies of the free space diagram to create 
a double free space diagram and searching for a monotone path which covers 
every point in P exactly once see Figure [TJ This algorithm can be used to show 
whether Sf(9P, dQ) < e and find the particular mapping(s) between dP and 
dQ. In turn this defines a placement of the diagonals, i.e., a mapping of the 
endpoints of the diagonals to endpoints of the corresponding image curves in Q. 

Subproblem (2) is solved by only considering paths through the free space 
diagram that map a diagonal d onto an image curve / such that /) < e. 

Naturally the particular placement of the diagonals determined in subproblem 
(1) could affect whether this is true. Therefore, they must check this for many 
paths in the free space diagram. Fortunately they can show that it is sufficient 
to only consider mapping a diagonal to an image curve which is the shortest 
path between the end points determined by the placement. 

Solving these subproblems generates a mapping between P and Q for s. This 
mapping might not be a homeomorphism but the authors show by making very 
small perturbations of image curves the mapping can be made into one. These 
perturbations can be arbitrarily small. Thus, because the Frechet distance is the 
infimum of all homeomorphisms on the surfaces, the Frechet distance is e. For 
simplicity we will refer to these generated mappings as homeomorphisms. By 
performing a binary search on a set of critical values they can use the above 
algorithm for the decision problem to compute the Frechet distance of P and Q. 

2.2 Shortest Path Edge Sequences 

Our algorithm extends the simple polygons algorithm to one for folded polygons. 
The idea of the algorithm is to subdivide one surface, P, into convex regions and 
pair those with corresponding regions in the other surface, Q. The difference is 
that those regions of Q are now folded polygons rather than just simple polygons. 
The authors of |Buc06j show that the Frechet distance of a convex polygon and 
a simple polygon is just the Frechet distance of their boundaries. Using almost 
the same argument we prove that it also holds for folded polygons, see Appendix 
E]for the full proof 




Fig. 1. The white areas are those in free space diagram. The surfaces are within 
Frechet distance e since there is a monotone path starting at the bottom of the free 
space diagram and ending at the top which maps every point on the boundary of P 
exactly once. This figure was generated using an ipelet created by Giinter Rote. 

Lemma 1. The Frechet distance between a folded polygon P and a convex poly- 
gon Q is the same as that between their boundary curves. 

As mentioned in Section 12.11 for the simple polygon algorithm it suffices to 
map diagonals onto shortest paths between two points on dQ. By contrast, there 
are folded polygons where a homeomorphism between the surfaces does not exist 
when diagonals are mapped to shortest paths but does exist when the paths are 
not restricted, see Figure [5] The curve s± is the shortest path between the points 
a and b but the curve S2 has smaller Frechet distance to d than si has. We must 
therefore consider mapping the diagonals to more general paths. Fortunately, we 
can show that these more general paths still have some nice properties for folded 
polygons. 




(a) (b) 

Fig. 2. An example where the image curve S2 with the smallest Frechet distance to a 
diagonal d is a non-shortest path si in Q. (a) overhead view, (b) sideview. 

Lemma 2. Let Q be a folded polygon, u and v be points such that u,v £ dQ, 
E = {eo 5 ei, . . . , e s } be a sequence of edges in the convex subdivision of Q, and d 
be a line segment. Given e > we can find a curve f in Q that follows the edge 
sequence E from u to v such that Sp(d,f) < s, if such a curve exists, in O(s) 
time. 

Proof. We construct a series F — FS e (eo, d), FS e (e±, d), . . . , FS e (e s , d), of 2- 
dimensional free space diagrams. Any two edges et and e^+i are on the boundary 



of the same convex polygon in the convex subdivision of Q so we can assume 
without loss of generality that / consists of straight line segments between the 
edge intersections. This is similar to the shortcutting argument used to prove 
Lemma 3 in [Buc06 . Thus, we only need to check the points where / crosses an 
edge of Q. For Sp(d,f) < e to be true, the preimages of those crossing points 
must be monotone along d. Let FS e (e i: d) be the projection of FS e (ei,d) onto 
d. Let F' = FS' £ (ei,d),FS' e (e 2 ,d),...,FS' e (e s ,d). 

To verify that the preimage points on d can be chosen such that they are 
monotone, we check the intervals of F . Specifically, for i < j, the point on d 
mapped to must come before the one mapped to 6j . This can be checked by 
greedily scanning left to right and always choosing the smallest point on d which 
can be mapped to some edge. A search of this form takes O(s) time. □ 

The dual graph of the faces of Q is acyclic. This implies that there is a 
unique sequence of faces through which a shortest path from u to v, where 
u, v £ dQ, must pass. Necessarily, there must also be a unique edge sequence 
that the shortest path follows. We call such an edge sequence the shortest path 
edge sequence. 

Lemma 3. Let d be a diagonal. If there is a curve f C Q with <5_f(cZ, /) < e then 
there is a curve g C Q which follows the shortest path edge sequence such that 
o~F(d,g) < e. 

Proof. Let Ef and E g be the edge sequences of / and g respectively. By definition 
the dual graph of the faces of Q is acyclic, so E g must be a subsequence of Ef. 
Eg induces a sequence of free space intervals. If there is a monotone path in the 
free space interval sequence induced by Ef, we can cut out some intervals and 
have a monotone path in the free space for E g . □ 
From Lemma 131 we just need to consider paths that follow the shortest path 
edge sequence. We refer to paths that follow this edge sequence and consist of 
straight line segments between edges as Frechet shortest paths. In addition, ,s in 
Lemma [2] is bounded by the number of edges in along the shortest path edge 
sequence between u and v, and E will be the shortest path sequence. This implies 
the following theorem. 

Theorem 1. Let Q be a folded polygon, u and v be points such that u, v 6 dQ, 
and d be a line segment. Given e > 0, we can in 0(1) time find a curve f in Q 
from u to v such that <5f(oJ, /)<£«/ such a curve exists. 

Suppose we have a homeomorphism between dP and dQ. The endpoints of 
the image curves must appear on dQ in the same order as their respective di- 
agonal endpoints on dP. The homeomorphism also induces a direction on the 
diagonals in P and on the edges in Q. Specifically, we consider diagonals and 
edges to start at their first endpoint along dP or dQ respectively in a coun- 
terclockwise traversal of the boundaries. We denote by D e the set of diagonals 
whose associated shortest path edge sequences contain an edge e C Q. Observe 
that pairwise non-crossing image curves must intersect an edge e in the same 
order as their endpoints occur on dQ. We refer to this as the proper intersection 
order for an edge e. 



2.3 Diagonal Monotonicity Test and Untangleability 

We now define a test between two folded polygons P and Q which we call the 
diagonal monotonicity test. For a given e this test returns true if the following 
two things are true. First, Sp(dP,dQ) < e. Second, for every diagonal di in the 
convex subdivision of P, the corresponding Frechet shortest path fi in Q has 
Sp{di, ft) < e. We refer to the class of mappings of the folded polygons generated 
by this test as monotone diagonal mappings. This is similar to the test used by 
|Buc06j except ours uses Frechet shortest paths instead of the shortest paths. 

Unfortunately, because the image curves of the diagonals are no longer short- 
est paths, they may cross each other and we will no longer be able to generate a 
mapping between the folded polygons which is a homeomorphism. For an exam- 
ple of this see Appendix [B] Thus the diagonal monotonicity test might return 
true when in fact a homeomorphism does not exist. We must explicitly ensure 
that the image curves of all diagonals are non-crossing. In particular, we refer 
to a set of image curves F = {/i . . . fk} as untangleable for e if and only if there 
exists a set of image curves F = {f 1 . . . f k } where /, and f i have the same end 
points on dQ, 5p(di,f i ) < e, and the curves of F are pairwise non-crossing. 
A homeomorphism exists between the folded polygons for e if and only if there 
exists a monotone diagonal mapping whose image curves are untangleable for e. 
This follows from the same argument used in Buc06 for simple polygons. 

As shown in Theorem [T] computing Frechet shortest paths instead of shortest 
paths does not increase the asymptotic run time. To optimize this e we can 
perform a binary search on a set of critical values. As in |Buc06] . the number 
of critical values is 0(m 2 n + mn 2 ). The three types of critical values between a 
diagonal and its corresponding path through Q are very similar to those outlined 
in the simple polygons algorithm. So, by following the paradigm set forth by 
|Buc06] . we arrive at the following theorem: 

Theorem 2. The minimum e for which two folded polygons P and Q, pass the 
diagonal monotonicity test can be computed in time 0(kT matr i X muit('m-n-) log(mn)) 

3 Fixed-Parameter Tractable Algorithm 

In this section we outline an algorithm to decide for a fixed mapping between the 
boundaries of a pair of folded polygons whether the image curves induced from 
the mapping are untangleable. From this we create a fixed-parameter tractable 
algorithm for computing the Frechet distance between a pair of folded polygons. 

3.1 Untangleability Space 

Let e be an edge in Q which is crossed by the image curves of h diagonals, 
d\ , . . . , dh . We assume without loss of generality that the image curves of the 
diagonals cross e in proper intersection order if, for all 1 < i,j < h where i < j, 
the image curve of crosses e before the image curve of dj crosses e. Let the 
untangleability space U e contain all fc-tuples of points on the diagonals which can 



be mapped to crossing points on the edge e within distance e and such that the 
crossing points are in the proper intersection order along e. U e can be shown to 
be convex yielding the following theorem. 

Theorem 3. U e (di, . . . , dk) is convex. 

This theorem can be proven by linearly interpolating between points in U e . 
For the full proof see Appendix [Cj 

3.2 Fixed-Parameter Tractable Algorithm 

We assume the complexity of k and I, the convex subdivisions of P and Q 
are constant. Checking for the existence of a set of image curves which are 
untangleable can be done by using the untangleability spaces of the diagonals. 
Assume we are given some homeomorphism between dP and dQ from which we 
get a placement of the diagonals. We first choose an edge in Q to act as the 
root of the edge tree that corresponds to the dual graph of Q. We propagate 
constraints imposed by each untangleability space up the tree to the root node 
to determine if the set of image curves induced by the placement of the diagonals 
is untangleable. 

The untangleability space of an edge e, U e , contains exactly 
those sets of points on the diagonals in D e which can be mapped 



to the edge in the proper intersection order. The point chosen in 
U e imposes a constraint on what points may be chosen in other d] 
untangleability spaces. In particular the corresponding points on 
all of the diagonals must be monotone with respect to their edge 



sequence. We define C(U ei ) as the Minkowski sum of U ei with 

a ray in the opposite direction of the constraint on each of the Fig. 3. 

diagonals in D ei , see Figure [3] In Figure [3] U e is shown in white 

and C(U e ) is the union of the white and light gray portions. The 

direction of this constraint depends on which side of the edge the next edge 

is. C(U ei ) contains exactly those sets of points on the diagonals not excluded 

from having a monotone mapping with U ei . 

We define for every edge e a /c-dimensional propagation space P e . If e is a leaf 
in the tree, then P e = U e . Otherwise, define 

p e = u E nc(p ei )n...nc(p ej ) 

where e is the parent of the edges e\, ■ • ■ , ej. C(P e . ) contains only those points 
that are not excluded by the constraints of the tree rooted at ej from being used 
to untangle on the parent of ej. The propagation space for the root will be 
empty if and only if this set of image curves are not untangleable. From our 
assumptions, the propagation space of the root can be computed in constant 
time as the intersection of semi- algebraic sets [Bas06 . Let F(k, I) be the time 
complexity of computing this intersection. 

Consider two different mappings between dP and dQ. These determine differ- 
ent placements of the diagonals. If all of the image curves of all of the diagonals 
have the same shortest path edge sequence in both of the mappings the test 



will return the same result. Thus, we only need to test paths through the free 
space diagram which cross the diagonals and edges in a different order. The free 
space diagram for dP and dQ contains 2k vertical line segments that will each 
contribute 0(kl) different mappings of the diagonals and edges. Hence, there 
are 0((kl) 2k ) paths through the free space diagram which we need to test. For 
each of these we can check whether a global untangling exists as described above 
in constant time. Similar to the algorithm for polygonal curves |Alt95] we can 
perform Cole's [Col87 technique for parametric search |Meg83| to optimize the 
value of s. For convex subdivisions with constant number of edges I and diagonals 
k, this yields a fixed-parameter tractable algorithm with runtime polynomial in 
to and n. 

Theorem 4. We can compute the Frechet distance of two folded polygons in 
time 0((F(k,l)(kl) 2k + kT matrixmu u{mn)) log(mn)). 

4 Constant Factor Approximation Algorithm 

In this section we present an approximation algorithm based on our diagonal 
monotonicity test to avoid tangles altogether. First we demonstrate the following 
theorem: 

Theorem 5. // two folded polygons, P and Q, pass the diagonal monotonicity 
test for some e, then 5f(P, Q) < 9e. 

Proof. Consider the image curves of the diagonals of P found by performing our 
diagonal monotonicity test. To pass the diagonal monotonicity test there must 
be a homeomorphism between dP and dQ. The image curves of the diagonals 
will be mapped in the proper order along the boundary of Q. Therefore, if a pair 
of image curves cross in Q they must do so an even number of times. 

Take two consecutive points u, v on the boundary of P that are connected 
by a diagonal d = uv, and consider the convex ear of P that d "cuts off". Let 
d be the image curve of d in Q. In order to create a homeomorphism between 
P and Q, d should cut off an ear of Q which can be mapped to the ear of P. 
Unfortunately, some image curves may cross this d and cause tangles. Consider 
the arrangement of image curves in Q and let d be the highest level of this 
arrangement closest to the top of the ear (dQ), such that d connects the image 
points u and v on the boundary of Q. 

Observe that if an image curve d 1 crosses d from below then that intersection 
point a has a pre- image on both d and d\. These points a on d and a\ on d\ can 
be no more than 2e apart since they both map to the intersection within distance 
epsilon. In addition, d x must cross back below d eventually since all image curves 
which cross do so an even number of times (take the first such occurrence after 
the initial crossing). The preimage points b and b\ of this second intersection b 
are also no more than 2e apart. Since both d and d are line segments, every 
point on the line segment ab on d is < 2e distance from some point on the line 
segment ai&i on d\. For an approximation of 3e we can map a point on d to its 



Fig. 4. (a) the intersections between d and di imply that ab can be mapped to the 
region of di between a and b within Frechet distance e (b) this is the preimage 
of the intersection between di and d 2 on d, (c) an example of the preimages of two 
intersections occurring out of order of d 

corresponding point on d\ and then to where that point maps on d 1 , see Figure 
[4ja). Hence, the diagonal d can be mapped to an image curve on or above d 1 
within Frechet distance 3e. 

If this image curve d 1 then crosses another image curve d 2 this argument 
above cannot be just repeated because the approximation factor would depend 
linearly on the number of image curves which cross each other. The preimage 
points on d of such an intersection not involving d are separated by at most 
6e. This is because both diagonals involved in the intersection have a 3e corre- 
spondence between the region of them mapped above d and d. If the preimages 
are in order there is no problem. If they occur out of order they cause a mono- 
tonicity constraint. Fortunately, we can collapse this region on d to the leftmost 
preimage with 6s and then map it to the corresponding point on d 2 in 3e for a 
total of 9e, see Figure 0|b). If the preimage points in this example were reversed 
they would be in order. 

Thus, we can approximate away the monotonicity constraint of single inter- 
sections with 9e. We must also verify that if the preimages of single intersections 
occur out of order it does not effect our approximation, see Figure HJc). Due 
to lack of space the discussion of these technical cases has been moved to the 
appendix, see Appendix [Di 

From these we get that 5p(d,d ) < 9e. Now collapse the ear we initially 
selected in P to d. Likewise in Q collapse the corresponding ear to d . This is 
okay because d is above all of the other image curves in Q. This pairs the ear we 
cut off of P with the part of Q above d which is a folded polygon (it could be 
simpler than a folded polygon but we know it's no more complex than that). By 
Lemma Q] the ear of P and the folded polygon above d must be within Frechet 
distance 9e. 

Choose another ear in P. We can repeat the above arguments to remove this 
new ear and its corresponding ear in Q. The dual graph of P is a tree. Each 
time we repeat this argument we are removing a leaf from the tree. Eventually, 
the tree will contain only a single node which corresponds to some triangle in P 
which we map to the remainder of Q. □ 

As a direct consequence of Theorems [5] and [2] we get the following theorem: 



Theorem 6. We can compute a 9 -approximation of the Frechet distance of two 
folded polygons in time 0(kT matrixmu i t (mn) log(mn)). 

5 Axis-Parallel Folds and Distance 

The counter example from Appendix [Bl works for all L p except for L^. In this 
section we outline a special case where using the metric guarantees that if 
a pair of folded polygons pass the diagonal monotonicity test for e their Frechet 
distance is no more than e. Specifically, if all of the line segments in the convex 
subdivision of the surfaces are parallel to the x-axis, y-axis, or z-axis, we show 
that it is sufficient to use shortest paths instead of Frechet shortest paths. Since 
shortest paths never cross we can use the simple polygons algorithm in [Buc06 
to compute the Frechet distance of the surfaces. We first prove the following 
lemma. 

Lemma 4. Let R be a half-space such that the plane bounding it, OR, is parallel 
to the xy-plane, yz-plane, or xz-plane. Given a folded polygon Q with edges par- 
allel or perpendicular to the x-axis and points a,b £ Q (1 R, let f be a path in Q, 
which follows the shortest path edge sequence between a and b. If P is completely 
inside of R so is the shortest path f between a and b. 

For the lemma to be false there must exist a Q, R, and / which serve as a 
counter example. There must be at least one edge ej in Q such that / n ej E R 
and / n ej $ R. In particular, let ej be the first edge where this occurs along 
the shortest path edge sequence. First consider a Q where all of the edges of it 
are perpendicular to OR. A line segment in the shortest path / connects the 
endpoints of two edges in Q. Let a and ek be the edges that define the line 
segment in / that passes through ej. We now consider several cases in how 
those edges are positioned. 

Case (I) occurs when ek is completely outside of R, see Figure [5ja). While 
this does force / to cross ej outside of 7?, there is no / which can pass through ek 
while remaining inside of R. Because Q is a folded polygon any path between a 
and b must path through the edges in the shortest path edge sequence including 
ek- Thus no / can exist entirely within R. 

Case (II) occurs when part of is in R and / crosses it in the part in R, 
see Figure [5fb). In this case / does not cross ej outside of R. 

Due to space limitations the discussion of the remaining cases has been moved 
to the appendix, see Appendix [E] Each of the remaining cases can be reduced 
to these first two. Using this lemma we can prove the following theorem: 

Theorem 7. The Frechet distance between two surfaces, both with only diag- 
onals/edges parallel to the x-axis, y-axis, or z-axis, can be computed in time 
0(kT matrixmult (mn) log(mn)). 

Proof. Let Q be a folded polygon, d be a diagonal, and / be the shortest path 
between points a and c on dQ. Using Lemma U we prove that if there exists a 
Frechet shortest path / between points a and c such that Sp(d, f) < e, then 
S F (d,f')<e. 




Fig. 5. (a), (b) are examples of case (I) and case (II). (c) example intervals for the 
two different paths, (d) together the edges ei and ei cause a monotonicity constraint. 

Minkowski Sum Constraints Since we are using the distance, the unit 
ball is a cube with sides of length 1. The Minkowski sum of a diagonal d in P 
and a cube of side length e yields a box. Points in the diagonal d can only map 
to points in this region. It can be defined by the intersection of 6 half-spaces; 
all of these have boundaries parallel to either the xy-axis, the xz-axis, or the 
yz-axis. Thus, from Lemma^l we know that if any path through Q is completely 
within this box, then the shortest path / will be, too. This means that for each 
edge ei on the shortest path edge sequence / D &i is within distance e of some 
non-empty interval of d. 

Monotonicity Constraints For the shortest path / between the boundary 
points to have Sp(d, / ) > e, at least two of these intervals must be disjoint and 
occur out of order along d, see Figure [5jc). Such a case introduces a monotonicity 
constraint on e. If no such intervals existed then we could choose a monotone 
sequence of points along d such that each point is within distance e of an edge 
and the sequence of edges they map to would have the same order as the shortest 
path edge sequence showing that 8p{d,f ) < e. 

Let e\ and e2 be two edges along the shortest path edge sequence for which 
such bad intervals occur. Let p\ and P2 be points on the shortest path where it 
intersect edges e\ and e2 respectively. Let q\ and q2 be the same for /. Finally, 
let d r contain all of the points on d which are within distance e of the point r. 
Since Sp(d, f) < e, d qi and d q2 must overlap or occur in order along d. 

Let i?i be the half-space whose bounding plane contains qi and is perpen- 
dicular to d. likewise let i?2 be the half-space whose bounding plane contains q^ 
and is perpendicular to d, see Figure [5jd) . Let R\ extend to the left along d and 
i?2 extend to the right along d. 9i?2 must occur before dR\ along d or the edges 
are in order and no monotonicity constraint is imposed. Assume R\ encloses all 
of / between a and &2- If it does not we can choose a new edge between a and 
e 2 to use as e\ for which this is true. Doing so only increases the monotonicity 
constraint. Likewise we can assume R2 encloses all of / between e2 and c. 

Assume a, b, and c lie on / . Specifically, let a and c be the end points of / on 
dQ. Naturally, a shortest path must exist between a and c and it must contain 
at least one point in R\ n R% which we call b. f follows the shortest path edge 
sequence between a and c, so it must also cross all of the edges in the shortest 



path edge sequence between a and b. Therefore, to show that P2 is inside of R\ 
we can directly apply Lemma |4] to the points a and b. A similar method can be 
used for ei with points b and c to show p\ is inside R2 ■ Since d qi and d q2 overlap 
or are in order, d Pl and d P2 must as well. Therefore, 5p(d, f ) < e and shortest 
paths can be used for this variant of folded polygons instead of Frechet shortest 
paths. Because we are using shortest paths we can just use the simple polygons 
algorithm. This yields Theorem [71 □ 

6 Future Work 

The constant factor approximation outlined in Section [4] can likely still be im- 
proved. Specifically, we consider only the worst case for each of the out-of-order 
mappings which may not be geometrically possible to realize. In addition, we 
currently approximate the Frechet distance by mapping image curves one-by-one 
to the top of the arrangement of other image curves. It would of course be more 
efficient to untangle image curves by mapping them to some middle curve rather 
than forcing one to map completely above the others. 

Finally, while the problem of untangling seems hard, it is also possible that a 
polynomial-time exact algorithm could exist. The acyclic nature of our surfaces 
seems to limit the complexity of our mappings. The methods used to prove that 
computing the Frechet distance between certain classes of surfaces is NP-hard 
in [BuclOj are not easy to apply to folded polygons. 
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A Proof of Lemma Q] 



Proof. Assume that the Frechet distance between the boundary curves is e. 
Now consider mapping the diagonals in the convex decomposition of P to Q. 
Because Q is convex, we can connect the images of the endpoints of each of 
the diagonals with a line segment. These line segments subdivide Q into a set 
of convex polygons each of which are paired with a convex polygon in P. The 
boundaries of the paired polygons must be within Frechet distance e of each 
other, consequently, by [Buc06] . the paired convex polygons must be in Frechet 
distance e. 

Combining all of these mappings, we have a homeomorphism between P and 
Q. Thus the Frechet distance of P and Q must also be < e. □ 

B Counter Example 

In this section we show a counter example to the claim that it is sufficient that 
a pair of folded polygons pass the diagonal monotonicity test for e to guarantee 
that they are within in Frechet distance e. Let P and Q be the surfaces shown 
in Figures [^a) and^b) with e = 1. 




Fig. 6. (a),(b) The folded polygons P and Q, (c) P and Q overlaid such that the zigzag 
portions of the surfaces outline a cube at their center, (d) The cube has sides of length 
one. The marked points must be used in the mappings for e = 1 but they can be shown 
to be mutually exclusive 

The surfaces are the same except rotated by 90 degrees. Each surface consists 
of two parallel 2 by 3 rectangles which are connected together by a smaller slanted 
rectangle. One of the large rectangles is at height 1 while the other is at height 
0. Figure WLc) shows how the two surfaces are overlaid. 

We can prove that the boundary curves have Frechet distance 1 to each other. 
Pair up each side of P with one in Q as indicated in Figure [T^b) . These paired 
portions of the boundary have Frechet distance 1 to each other. The case of 
mapping A to A is given in Figure Efc). For the beginning and ending of the 
sides the curves either overlap or are exactly distance 1 apart with one directly 
above the other, so the mapping used is obvious. Matching the flat middle part 
on A to the zigzag in A is a bit harder. The key insight is to observe that the 
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Fig. 7. (a) Sideview of the surfaces, (b) Top view of P and Q showing which sides are 
mapped to which, (c) A side of P is mapped onto a side of Q. ai and as well as a 2 
and o 2 occupy the same points in space and a separation is shown for clarity. 



point a 2 is within distance 1 to 02, 03, and 04. The entire zigzag can be mapped 
to this point. The portion after the point a 4 can then be mapped straight down. 

Consider the intersection between P and Q. At the center of the folded 
polygons is a cube with side length 1, see Figures [5{c) and[SJd). The point p\ 
on d\ can map both to the point p 2 on edge e 2 and to the point r on edge e\ 
within distance 1. Thus, there exists an image curve /1 such that Sf(cIi, f\) < 1. 
Likewise, the point g 2 on d 2 can map both to the point s on edge e 2 and to the 
point gi on edge e\ within distance 1, so there exists an image curve f% such 
that (5f(^2, / 2 ) < 1. Therefore, P and Q pass the diagonal monotonicity test for 
1. 

Notice that p\ is the only point on d\ which can map to e%. Given the direction 
of di, fi starts at the top of Q and must cross edge e 2 then e%. Thus no points 
after p\ on d can mapped to e 2 so pi is the only reachable point on e 2 . Likewise, 
s is the only point on e 2 the diagonal d 2 can map to. s comes before p 2 along 
e 2 so the image curve of <i 2 must cross e 2 before the image curve of d%. This 
is different from the proper intersection order that comes from the boundary. 
Therefore, even though P and Q passed the diagonal monotonicity test for 1, 
Sf(P,Q) > 1. 

C Convexity of the Untangleability Space 

Proof. To construct a proof we first define U e more formally. 

U e (di, . . . ,d h ) = {a = (ax,...,a h )\(3pi,...,Ph G e|(Vi,j|(a 4 G di, 1 < i < j < h) -> 

(IK - Pill < e,Pi <Pj)))} 

where pi < pj means that the point pi is no further than the point pj along 
the edge e. If U e is not convex then there are points a = {a\, . . . , a/,), c = 
(ci, . . . , Ch) G U e such that there exists a point b — (b\, ... , bh) on the line 
segment between a and c where b U e , see Figure 15^. The point & is on the 
line segment between a and c so there exists an x between and 1 such that 
b = xa + (1 — x)c. Likewise for each point bi £ 6, bi — xai + (1 — x)ci. Let Pi 
and Ti be the points on edge e mapped to by pointa <Zj and respectively. Let 



q. L = xpi + (1 — x)ri. Since a,c € U e , ||<Zj — Pi\\ < £ and ||cj — T"» [| < e. It follows 
from these that — < e, see Figure [5Jd. 

||<ii < £ z||a» -Pi|| < 
=>■ ||xai — xp, || < 

||ci - nil < e ^ (1 - x)||cj - ri|| < (1 - x)e 

- x) Ci - {1 - x)n\\ < {l-x)e 
=>■ ||(1 — x)ci ~ (1 — x)ri||||xaj — xpi\\ < (1 — x)e + xe 
=>■ \\(xai + (1 — x)ci) — (xpi + (1 — £c)rt) || < e ~ xe + xe 
=> \\bi -qi\\ <£ 

Since a,c E U e , Pi < and r^ < rj. It follows from a proof similar to the 
previous one that % < qj, see Figure |5b. 

Pi < P 3 iCPi < £Pj 

< rj =>■ r. ( < rj 

=^ (1 — a;)r,; < (1 — x)rj 

=> xpi + (1 — x)rj < xpj + (1 — x)rj 

=> q i < q 3 

Thus, for all where < i < j < h, ||6j — <Zi|| < £ and q,; < qj. From this 
it follows by definition that b s U e . Therefore U e is convex. □ 




Fig. 8. (a) a,c £ [7 e but b ^ U e - (b) both and Ci being within distance e of their 
respective points on e necessitates that 6; is within distance e of its own. (c) The 
crossing points of diagonals di and dj being in order along e for both points a and c 
necessitates that they are also in order for b 



D Constant Factor Approximation Order of Regions 



We now need to show that the points that the intersections are mapped to do 
not occur out of order. This could happen when the preimage regions on d occur 
out of order, see Figure IHc). In this case, collapsing to the leftmost point of a 
region for a single intersection would not be sufficient. The two collapsed points 
would still be out of order along d. 

We consider two cases. For each of these cases we assume that the intersection 
points a and b are at on the top level of the arrangement of the image curves. 




Fig. 9. (a) If the image curves intersect before crossing d , the preimages of both 
intersections must completely occur on a region of d smaller than 6e, (b) If the image 
curves intersect d before each other, the preimages of the intersections must occur in 
the proper order 



Case I 

Sec figure Figure IHJa). Consider the image curve of the first intersection which 
is sloping toward d (d 1 in example) and the image curve of the second intersec- 
tion which is sloping away from d (d 4 in example) intersect each other before 
intersecting d (moving left-to-right). Assume d 1 first intersects d 2 and then <i 4 , 
and assume d 4 first intersects d 2 and then d 3 . 

In this case we still have a correspondence between the two diagonals and 
d (they are still above d and thus a mapping exists). The preimages of the 
intersection points on d are restricted to being outside the regions from Figure 
EJc). The region they enclose must be < 6e in width. Thus we can just map the 
whole region to the leftmost point. 

Case II 

See Figure [SJb). The image curves cross d before they intersect each other. The 
order the curves intersect in places a restriction on how far back the preimage 
of intersection 2 can be. Specifically, in this case, the preimages of the two 
intersections must occur in order on d so there is no new restriction imposed by 
them. 



Induced Mapping 



One of these two cases must occur. In particular, the case where the image curves 
cross earlier or later than specified in Case I cannot happen. In this hypothetical 
case one of the image curves would be forced above the other one. Assume d 1 is 
the one forced above. In that case d 1 would need to pass above the intersection 
b in order to not cross d 4 . We assumed earlier that b was at the highest level 
of the arrangement so this is a contradiction. The same argument works if d 4 is 
the image curve forced above. 

From these two cases we can see that monotonicity constraints will not occur 
between different intersection preimages on d. Thus, we can collapse each of 
them to a single point using the 9e approximation described above without the 
collapsed points occurring out of order. 

E Remaining Cases for Lemma 31 

Case (III) occurs when part of e k is in R and / crosses it in the part outside of 
R, see Figure [TOT a). The path / can cross the part of e% inside R. This case can 
be repeated many times but eventually, / ends at b so it needs to cross back 
into R. To do this it must cross an edge on the opposite end point. If that end 
point is in R the entire path can be shortcut similar to case (II) and if the point 
is not in R then as in case (I) there does not exist an / which can cross e/ and 
remain in R see Figure fTOT b') . Thus case (III) also leads to a contradiction. 




Fig. 10. (a) the edge e k +i forces / out of R on e 3 -, but for this to happen e k +i must 
be completely outside of R, thus no / can cross e k +i- (b) / is not forced out of R on ej 
but rather can follow the path / . (c) the path can continue to curve away from R. (d) 
in this case, while the shortest path goes out of R eventually there must be an edge e; 
where it reverses direction since / ends at b which is in R. Such an edge cooresponds 
to case (I) . 

Next we consider the two cases that arise from adding a perpendicular fold 
e k+ \. Specifically, we will consider where the first perpendicular fold after ej is 



placed. Such a fold runs parallel with dR so it must be either entirely inside R 
or entirely outside. 

Case (IV) occurs if the edge ek+i is outside of R. In this case there does not 
exist an / which can pass through eu+i and be entirely in R, see Figure [TUT c). 

Case (V) occurs if the fold et+i is inside of R. In this case the shortest path / 
will no longer be forced outside of R on the edge ej +1 . Instead, it goes through R 
to the edge efc+i, see Figure HUTd). The only way that / could not pass through 
R would be if one of the edges between ej and were entirely outside R and 
that is ruled out by case (I) 

Thus no such Q, R, and / can exist. 



